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Abstract 

Using the CLEO II detector at CESR, we measure the r/' energy spectra in T(15) decays, that 
we compare with models of the g'g*g form-factor. This form-factor especially at large g' energies 
may provide an explanation of the large rate for B Xgg'■ Our data do not support a large 
anomalous coupling at higher cp' and thus the large g' rate remains a mystery, possibly requiring 
a non-Standard Model explanation. 
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I. INTRODUCTION 


There are several interesting, nnexplained phenomena in B decays. First of all, the total 
prodnction of charm and charmoninm seems about 10% low especially when coupled 
with a B semileptonic branching ratio of (10.4±0.3)% j^- Secondly, CLEO observed a very 
large rate of rj' in the momentum range from 2 to 2.7 GeV/c with a branching fraction of 
(6.2 ± 1.6 ± l.Sl?;®) X 10“^ j^. The BABAR experiment has conhrmed this large rate 0. 
The production of r]' mesons is believed to occur dominantly via the b ^ sg mechanism, as 
strongly suggested by observation of the two-body decay B rj'K. One explanation of the 
large g' rate is that the b ^ sg rate is not 1% as expected in the Standard Model, but is 
enhanced by new physics to be at the 10% level. This would also explain the charm deficit 
problem. 

An alternative explanation is that of an anomalously strong coupling between the rj' and 
two gluons iiQ- The process b ^ sg followed by the two gluon coupling to the rj' is 
shown in Fig. ^ 



FIG. 1: Diagram for b —> sgg'. 


Experimentally, the hadronic mass associated with Xg sometimes is a K, ~10%, and even 
more rarely a K*, ~1%; in fact, most of the rate has the mass of the Xg system larger than 
1.8 GeV. 

Since the g' is mostly the flavor singlet gi, as the g — g' mixing angle is between 10° — 20°, 
the effective g'g*g coupling can be written as 

H{q‘^)eai3^uq°‘k^£ie2, (1) 

where q = pb — pg is the four-momentum of the virtual hard gluon {g*), k is the four- 
momentum of the soft “on-shell” gluon (g), and H{q^) is the g*gg' transition form factor. 

Ghen and Kagan [sj have shown that the region of the q^ relevant in the process b —> sgg' 
can also be accessed in high energy g' production in T(IS') decay. Thus constraints can be 
put on the H{q^) from the g' spectrum in T(IS') —ggg decays. H{0) is found from the rate 
of JI'll) —> -gg' decays as ~1.8 GeV“^. 

Three choices for the form factor shape H (g^) are shown in Fig.|21 (a) a slowly falling form 
factor from Hon and Tseng H{q^) = 2.1 GeV ^ag{q^)/ag{ni?^,)] (b) a rapidly falling form 
factor representative of perturbative QGD calculations, H{q^) = 1.7 GeV“^m^,/(g^ — rri^,) 
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at > 1 GeV^; (c) an intermediate example with H{q^) cx l/(g^ + 2.2^GeV^) [^. In (b) 
and (c) the form factor at ~ has been matched onto the valne given in (a), which is 
hxed by the QGD anomaly j^. The parametrization of the form factor in (b) follows from 
a simple model in which the rj' is conpled pertnrbatively to two glnons through quark loops 
[^. With the choice H{0) = 1.7 GeV“^ it compares well with the perturbative QGD form 
factors obtained by other authors 
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FIG. 2: Three choices for the form factor plotted against (a) the slowly falling form 

factor, (b) a rapidly falling form factor representative of perturbative QCD calculations, (c) an 
intermediate example (adapted from j^). 


We will compare the theoretical predictions for H{q^) with data taken on the T(IS') 
resonance with the GLEO II detector at the GESR storare ring. Some information on this 
topic has been extracted by Kagan from ARGUS data 


II. DATA SAMPLE AND ANALYSIS METHOD 

In this study we use 80 pb“^ of GLEO II data recorded at the T(IS') resonance (9.46 
GeV), containing 1.862 x 10® T(IS') events. We also use off-resonance continuum data 
collected below the T(4S') resonance (10.52 GeV) with a total integrated luminosity of 1193 
pb“L 

The theoretical predictions referred to in this paper are made for T(IS') decays into three 
gluons {ggg). In order to compare our measurement to them we have to correct for the 
T{1S) ^ —>■ qq contribution, whose size is given by 

i3(T(l^) ^ 7* ^ gg) = R . i3(T(l^) ^ = (8-8 ± 0.3)%, (2) 

where = 3.56 ± 0.07 [l^ and the i3(T(lS') —/i■*■/!“) is taken as (2.48 ± 0.06)% [^. 

Although several processes can contribute to inclusive rj' production in T(IS') decays, 
it is believed that the soft processes including fragmentation populate only the low g^ or 
equivalently the low Z region, where 

Z = En'/Ebeam = ‘2:Erj'/M{T{lS)). (3) 
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Thus in the large Z region significant rj' production would indicate a large ri'g*g coupling. 

The CLEO II detector, described in detail elsewhere Q had a high resolution elec¬ 
tromagnetic calorimeter comprised of 7800 Csl crystals surrounding a precision tracking 
system. 

We detect g' mesons using the decay channel: g' with a branching fraction of 

44%, and > 77 with a branching fraction of 39%. We identify single photons based on 
their shower shape and the non-proximity of charged tracks. Those photon pairs within the 
“good barrel” region of the detector, | cos^l < 0.707 (where 6 is the angle with respect to the 
beam), that have invariant masses consistent with the g mass within 3 standard deviations 
are constrained to have the invariant mass of the g. For g mesons coming from low energy g' 
candidates (Z < 0.5) the background from 71° decay is large, and thus the candidate photons 
are also required not to be from a possible 7r° decay. We then add two opposite sign pions 
and form the invariant mass. 

The g7T~^7r~ invariant mass spectra are shown in Fig. |21for T(IS') and for off-resonance 
continuum data. The spectra are fit with a Gaussian function for signal and second order 
polynomial function for background. The numbers of reconstructed g' are extracted from 
the fit. We find 1486 ± 137 g' from the T(IS') data, and 4062 ± 174 g' from the off-resonance 
data. 
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FIG. 3: The g7r~^7r~ invariant mass spectrum reconstructed from T(15) data (left), and off reso¬ 
nance data at 10.52 GeV (right) fit with Gaussian functions for signal and second order polynomials 
for background. 

To measure the energy spectrum we reconstruct g' candidates in Z intervals. We choose 
the Z steps as 0.1. The invariant mass spectra are fit with the same functional form as 
used for Fig. 01 Here we fix the mass of the g' to our average value over all Z; the Monte 
Carlo simulation shows that the mass measurement should be independent of g' energy. We 
extract the width of the signal Gaussian distribution from Monte Carlo simulation for each 
Z bin and perform a smooth fit as a function of Z. The smoothed values are used in the fit 
as fixed parameters. The Z dependent mass spectra are shown in Fig. (Hand Fig. 01 

for T(IS') and off-resonance data, respectively. 

In order to extract decay rates we need to correct our raw event yields by efficiencies. 
These may not be equal for different intermediate states, i.e. qq versus ggg. The hadronic 
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FIG. 4: The r]7r~^7r~ invariant mass spectra in different Z ranges reconstructed from T(IS') data, 
fit with a Gaussian function for signal and a second order polynomial for background. 

events at T(IS') energy arise from different sources: about 4 nb is qq from continuum e''“e“ 
collisions, about 2 nb from T(IS') —*• 7 * —>• qq, 18 nb from ggg, and 0.5 nb from jgg from 
the T(IS'). The hrst two have same event topology and reconstruction efficiencies. We use 
the qq Monte Carlo generator to simulate these events. The 'jgg events are similar to that 
of ggg and have a relatively small cross-section; thus we treat them the same way as ggg 
events. We use the ggg Monte Carlo generator to simulate this part. 

We rely on off-resonance continuum data to estimate the qq contribution in T(IS') data. 
However, the continuum data were taken for continuum subtraction in T(4S') studies. The 
center of mass (CM) energy (10.52 GeV) is close to T(4S') mass (10.58 GeV), but more than 
1 GeV higher than T(IS') mass (9.46 GeV). The difference of reconstruction efficiency due 
to this energy difference is not negligible. We thus use different qq simulations for continuum 
data and T(IS') data. 

The energy difference also affects the Z spectrum of g' from continuum Monte Garlo as 
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FIG. 5: The rjir'^TT invariant mass spectra in different Z ranges reconstructed from off-resonance 
data, fit with a Gaussian function for signal and a second order polynomial for background. 


shown in Fig. EKa). The solid line is the E^//Ebeam distribntion for the T(IS') data (9.46 
GeV) and dashed line for the continunm data (10.52 GeV). The low limits are 0.202 and 
0.182 respectively. The discrepancy is signihcant, especially at low energy. In order to use 
our continuum data at 10.52 GeV we need to map it to 9.46 GeV. To do so we rely on the 
continuum Monte Garlo. We take the two Monte Garlo rj' shape distributions at 10.52 and 
9.46 GeV, denoted by '^’ 10 . 52 ( 2 :) and 'P 9 . 46 (^) and numerically integrate them to satisfy the 
relation: 

Vio.Mdz = V,A,{z)dz , (4) 

Jo Jo 

where Z(q 52 is fixed and a value for is determined. The data points on Fig. EKb) show 
the difference in Zg 45 —Z'^q 52 as a function of Z'^q 52 (or equivalently Zq in following function). 
We fit the points with a fourth order polynomial function to dehne the mapping analytically 
as 

Z = -0.215 X 10-2 + 1.2238 Zq - 0.6879 Z^ + 0.8277 Zj) - 0.3606 Z^ . (5) 


7 




























































The simplest mapping would be a linear conversion Z = 0.025 + 0.975 x Zq, shown as dotted 
line in Fig. EKb). We use this alternative to estimate the systematic uncertainty due to the 
mapping. 

That this mapping works is demonstrated in Fig. |Ul-(a), where the spectra shown as open 
circles is the mapped spectrum according to Equation El It overlaps well with the Monte 
Carlo spectrum generated at 9.46 GeV. 



FIG. 6; a) The Z = E^t/Eheam distributions from Monte Carlo simulation. The solid line is 
the Z=Erii /Eheam spectrum for an energy of 9.46 GeV, the dashed line is the spectrum for 10.52 
GeV and the open circles are the mapped spectrum from 10.52 GeV. b) The data points show the 
difference in the Z values at 9.46 and 10.52 GeV as a function of the Z value at 10.52 GeV. The 
solid curve is a fit to a fourth order polynomial. The dotted line shows the mapping of the linear 
conversion. 

The 7]' production rate is smaller at 9.46 GeV because of less available energy. From the 
qq generator we found that the production rate is 93.6% that of 10.52 GeV. This factor is 
also considered in estimation of the r]' production from qq events. 

The mapping for continuum data is derived from the model-dependent Monte Garlo 
spectrum. If the real data and the Monte Garlo are very different then the systematic 
uncertainty due to this mapping could be large. To check this, we compared the measured 
E,,//Ebeam spectrum with the generated spectrum. Fortunately, the spectra agree reasonably 
well and the systematic uncertainty due to this source is negligible. 

We now turn to estimating the detection efficiencies. Shown in Fig. [7| are the efficiencies 
estimated with different models and different energies for a) without the n° veto, and b) with 
the 71° veto. In the real data we applied 7r° veto to rj' candidates with Z < 0.5. Gomparing 
with the efficiency from 9.46 GeV qq events, the efficiency from ggg events is roughly 15% 
higher, and the efficiency from 10.52 GeV qq events is roughly 7% lower. The main source 
of such difference is the event shape. The ggg events are more spherical while the higher 
energy qq events are more jetty. 
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FIG. 7; The rj' reconstruction efficiencies as function of Z for different MC samples a) without a 
7 r° veto, and b) with a 7 r° veto in the photon selection. 

III. EXTRACTION OF THE r/' SPECTRUM FROM T(15) DECAYS 

The T(IS') data sample can be broken down into three parts as described in the previons 
section: 

Nall = Ny(^is)^ggg + iVT(lS)^gq + ^e+e-^qg • 

The hrst one has different reconstrnction efficiencies from the other two. For the contribntion 
from continnnm (e^e“ —> qq) events, we mnltiply the nnmber from off-resonance events at 
10.52 GeV, mapped nsing EqnationEl by a factor fa+e-^qq defined as: 

A(p )e-he-^qg-(9.46 GeV) = iV(V)e-He-^qq(10.52 GeV) X fa+e-^qq, (6) 


where 


fe+e-- 


80.4 1/9.46^ 


►<?<? 


X 

1193 

= 0.078 X 


1/10.522 

^ 9.46 

^ 10.52 


X 0.9356 X 


^ 9.46 

^ 10.52 


(7) 


where the first factor is the relative Inminosities, the second the energy sqnared dependence 
of the cross section, the third the relative r]' yield and e is the Z-dependent reconstrnction 
efficiency for qq events as shown Fig. [3 

We also want to evalnate the yield from T(IS') —> 7 * — qq. Since we know that 


= 0.555 ± 0.022 nb and (9.46 GeV) = 1.12 nb jl^, we derive the 


factor to be nsed in the Nr{is)- 


,qq estimation as: 


fr(lS)^qq — fe+e-- 


*qq 


X 




R ■ CTp 




= 0.0387 X 


^ 9.46 
^ 10.52 


( 8 ) 
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In Table HI we list the number of reconstructed r]' over all Z and in the high Z region for 
various T(IS') and continuum yields (only statistical errors are shown). Note that the total 
numbers of signal from T(IS') data and off-resonance data in this table are the sum of all 
Z bins derived bin per bin, as we need to use Z-dependent efficiencies. 


Sample 


All Z 

Z > 0.7 

T(IS') data 

1494 ± 120 

46.0 ± 8.1 

off-resonance 

4294 ± 130 257.1 ± 17.3 

T(15) - 

ggg 

972 ± 120 

13.9 ± 8.1 

T(15) - 

QQ 

173 ± 5 

10.6 ± 0.7 

Continuum qq 

349 ± 11 

21.5 ± 1.4 

T(1F) - 

ggg,QQ ii45 ± 120 

24.5 ± 8.1 


TABLE I: Number of reconstructed rj' from T(IS') and off-resonance data and the breakdown 
categories of T(15) data. Also listed are for samples with Z > 0.7. 

The measured T(IS') — rj'X branching fractions are listed in Table HTl both for Z > 
0.7 and for all Z. In the large Z region for 3 gluon decays, we do not have a statistically 
signihcant signal and thus derive a 90% conhdence level upper limit of B(T{1S) —>• ggg — >• 
rj' X) z>o.7 / B(T(IS) —>• ggg) < 3.4 x lO”"^. We describe the systematic errors below. 


Mode 




All Z 

Z > 0.7 

B{T{1S) 

^g'X) 



(2.8 ±0.4 ±0.2) % 

(3.1 ±0.9 ±0.3) x 10-^ 

B(T(15) 

ggg - 

g'X)/B{T{lS) - 

ggg) 

(2.8 ±0.5 ±0.2)% 

(1.9 ± 1.1 ±0.2) x 10“^ 

B(T(15) 

qq^ 

g'X)/B{T{lS) - 

> qq) 

(4.2 ± 0.2 ± 0.4) % (16.8 ± 1.1 ± 1.7) x lO”"^ 


TABLE II: Branching fractions of T(IS') to g' mesons, for all decays, three gluon decays and 
quark-antiquark decays for the entire g' energy spectrum and for Z > 0.7. The errors after the 
values give the statistical and systematic uncertainties, respectively. 


The sources of systematic uncertainties are listed in Table imi along with estimates of their 
sizes. The total systematic errors on branching ratios are ±10% for qq sample (independent 
of Z), ±11% for ggg sample at Z > 0.7, and ±8.6% for the rest. 

We also measure the differential branching fractions as a function of Z as shown in Fig. |S1 
In these plots only the statistical error is shown, which dominates the total error. 

We dehne three relevant differential branching ratio’s dn/dZ as: 


dnjggg) ^ djg(T(l^) ^ ggg g'X) 
dZ dZ X BiT{lS) ggg) 

dn{qq) _ dB{T{lS) ^ qq ^ g'X) 
dZ ~ dZ X B{T{1S) qq) 

dn{lS) dB{T{lS) g'X) 
dZ ~ dZ ■ 


(9) 


Listed in Table HU are the differential branching fractions in Z intervals for T(IS') decays 
to g' for ggg and qq subsamples and all decays. 
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dn/dZ 


Sources 

ggg {Z > 0.7) 

m 

All others 

Reconstruction efficiency of vr^ 

4.4 

4.4 

4.4 

Reconstruction efficiency of g 

5 

5 

5 

Number of g' from fit 

2 

2 

2 

Total number of T(IS') 

2.4 

2.4 

2.4 

B{g' —> 7r+7r“r/) 

3.4 

3.4 

3.4 

B(T(15) ^ 

- 

3.2 

- 

Ratio of integrated luminosity 

2.9 

1 

- 


3.6 

4 

- 

Z mapping 

6 

3 

3 

Total 

11 

10 

8.6 


t We use ^ {qq)) = (8.83 ± 0.28)%. 


TABLE III: The systematic uncertainties (in %) from different sources on the branching fraction 
measurements for the 3 gluon sample for Z > 0.7, the qq sample, and both the 3 gluon sample for 
all Z and the total T(IS') sample. 
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FIG. 8: The differential branching fraction dnIdZ as defined in context for a) T(IS') —> ggg —> g'X, 
b) T(1,S) -^qq^ g'X, and c) T(15) ^ g'X. 


In the Z spectrum of rj' mesons produced via ggg^ there is an excess above an apparent 
exponential decrease for 0.6 < Z < 0.7, corresponding to a recoil mass opposite the g' in the 
range 5.3 to 6.1 GeV jl^. However, a detailed study did not reveal any narrow structures. A 
possible explanation is that there is more than one process contributing to this distribution. 
We note also that the qq has much larger rates at high Z than ggg. 


IV. COMPARISON WITH THEORY AND CONCLUSIONS 

Fig. ini shows the Z spectrum of the g' measured in this paper compared with the spec¬ 
tra predicted by the three different models described above. The models are expected to 
dominate g' production only for Z> 0.7, with other fragmentation based processes being 
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z 

T(1S) ^ {ggg) T(15) ^ {qq) 

All T(15) 

0.2 - 0.3 

11164 ± 4471 

13205 ± 1253 

10503 ± 3740 

0.3 - 0.4 

11624 ± 1314 

11250 ± 685 

10716 ± 1099 

0.4 - 0.5 

3381 ± 558 

8898 ± 416 

3614 ± 467 

0.5 - 0.6 

1067 ± 300 

5030 ± 272 

1336 ± 251 

0.6 - 0.7 

963 ± 181 

2321 ± 166 

1011 ± 151 

0.7 - 0.8 

184 ± 92 

1116 ± 102 

252 ± 77 

0.8 - 0.9 

5 ± 50 

415 ± 59 

41 ± 42 

0.9 - 1.0 

31 ± 22 

153 ± 36 

40 ± 19 

0.7 - 1.0 

19 ± 11 

168 ± 11 

31 ± 9 

sum of all 

2842 ± 471 

4239 ± 153 

2751 ± 394 


TABLE IV: Differential branching fractions of r]' (xlO“^). The last two rows are total branching 
fractions. The branching fractions in columns 2 and 3 are normalized to the total branching fraction 
of T(15) —> {ggg) and T(IS') ^ {qq) respectively, while the last column is normalized to all T(IS') 
decay. The errors are statistical only, the systematic errors on the absolute normalization for 
column 1 is 8.6% for Z < 0.7, 11% for Z > 0.7, and 10% and 8.6% for columns 2 and 3, respectively. 


important at lower Z. The measurement strongly favors a rapidly falling dependence of 
the g*grj' form factor predicted by pQCD P, 10| , and ruling out other models. 

In conclusion, we have made the first measurement of the g' energy spectrum from 
T(IS') —> ggg decays. Our data are not consistent with an enhanced g'g*g coupling at 
large g' energies. Thus, the large observed g' yield near end point of the charmless B decay 
spectrum cannot be explained by a large g'g*g form-factor. Therefore, new physics has not 
been ruled out and may indeed be present in rare b decays. 
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z 


FIG. 9: The measured dn/dZ spectrum of T(15) ^ {ggg) g'X compared with theoretical pre- 
dictious. Showu iu dots are the measuremeut iu this study. Showu iu hues are differeut theoretical 
predictious: a) a slowly falliug form factor, b) a rapidly falliug form factor, aud c) iutermediate 
form factor |^. These predictious are valid ouly iu the regiou Z > 0.7. 
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